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We show that eigenstates |n) of the harmonic oscillator coupled to a linear Markovian bath demonstrate a 
non-trivial behavior in the dynamics of measures of non-classicality. Specifically, as the system undergoes 
decoherence, a time-dependent peak in non-classicality as a function of n emerges. We find this effect studying 
the dynamics of several non-classicality measures previously presented in the literature lOJ-Q] which compare 
quantum states to the set of all classical states. In studying these measures we introduce a novel set of classical 
states for the puipose of calculations which improve upon the results obtainable using these measures. In 
addition, following in the footsteps of [4], we show that the negative volume of the Wigner function agrees 
well with all the non-classicality measures while being otherwise calculationally significantly more tractable. 
Finally, we explore the dynamics of non-classicality of several other states. 

PACS numbers: 03.65.Yz 
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I. INTRODUCTION 



The macroscopic world is constructed of fundamentally quantum objects, but the observation of quantum effects requires 
either very small or very cold systems. That is, quantum effects are most easily visible in systems sufficiently isolated from their 
environments to prevent decoherence. It is generally more difficult to adequately protect larger systems from decoherence; this 
difficulty gives rise to the classical behavior of the macroscopic scale J5[]. Contemporary investigations of the transition from 
quantum to classical are motivated by fundamental considerations as well as practical issues of control and engineering. The 
transition is being investigated in nanomechanical systems JS, Q], mesoscopic systems such as Josephson junction devices |H], as 
well as cavity-QED (quantum electrodynamics) systems DSD, among others. This quantum-classical transition is controlled by 
internal system parameters and behaviors as well as (a) the relative size of h compared to the characteristic action of the system, 
(b) the strength of the system-environment interaction, and (c) the temperature and other characteristics of the environment, 
leading to a rich landscape of possible behaviors ifTolfTTll . 

The decohering harmonic oscillator has been the subject of many studies over the years lfT2l - [T6ll . While these studies are 
highly relevant for the fundamental understanding of quantum mechanics itself, it is not just a theoretical exercise since there 
have been proposals to encode qubits and qudits (relevant in quantum information theory) in the harmonic oscillator using finite 
superpositions of eigenstates [17]. In addition, experimental data now exists on the decay of these states where Brune et. al. 
successfully excited Fock states in a cavity-QED system, and monitored their decay using tomography II 1611 . 

In this paper we consider the quantum to classical transition in a decohering harmonic oscillator specifically studying the 
dynamics of measures of "non-classicality" in this system. There have been several careful studies of non-classicality measures 
in static, closed systems, and other detailed studies of decoherence in open systems; however, to the best of our knowledge there 
has been only one previous study of the dynamics of measures of non-classicality IU2I1 which looked at a restricted class of such 
measures on a restricted class of states. The study of these non-classicality measures in dynamical situations (that is, in open 
decohering systems) allows for deeper insight into the dynamics of decoherence while simultaneously providing a stronger test 
of the meaning and validity of these non-classicality measures. Our focus is on the harmonic oscillator eigenstates (the Fock or 
number states) and we also consider finite superpositions of these states. 

We show below that indeed these measures of non-classicality have non-trivial dynamics for an open system. Specifically, as 
the system evolves, a time-dependent peak in "non-classicality" as a function of eigenstate n emerges. In studying this behavior 
we also compare several measures of non-classicality and show that a somewhat intuitive measure (the negative volume of the 
Wigner function representation of the density matrix) works as well as or better than the other more formal measures. 

We start our presentation of these results in Section II with a new derivation of an analytical solution in the number basis of the 
behavior of the density matrix of a harmonic oscillator coupled to a linear Markovian bath. Although similar solutions have been 
previously derived, representing it in the number basis allows us to conveniently compute non-classicality as a function of time 
using several measures. Most of these measures are formulated in Hilbert-space, as opposed to characteristic or Wigner space. 
Thus the convenience comes from the fact that our solution allows for extremely fast calculation of the individual elements of 
the density matrix for any time and any bath temperature. In Section III we discuss in some detail how non-classicality has 
been defined and quantified. We first present different historical perspectives on the matter in order to gain intuition for the 
concept, and then go on to present four formal definitions of non-classicality from the literature which we believe provides a 
representative sample of the different techniques currently in use. 

Section IV is the heart of our discussion where we present results, starting with the static case; in the process of analyzing 
these, we are able to generalize previous results in that we find a novel classical basis that improves upon previous calculations 
for the static case (and is later of value for the dynamical case). We then present our results for decaying Fock states coupled 
to the vacuum, showing that during the transition from the quantum to the classical regime, a time dependent peak in non- 
classicality exists as a function of eigenstate n. We discuss how these results align with our intuition about classical-quantum 
correspondence. In addition, these results allow us to validate previous arguments that the negative volume of the Wigner 
function is a useful and meaningful measure of non-classicality, in that it agrees well with all the other non-classicality measures 
while being otherwise computationally significantly more tractable. The possibility that the Wigner function is a meaningful 
measure of non-classicality has current significance particularly because suggestions of how to directly measure the Wigner 
function, put forward two decades ago, are now being put into practice in many studies 118142 ill . This suggests that non- 
classicality could be easily measured in 'real time' without the need for detailed analysis. Indeed, since several of the measures 
of non-classicality we study require computationally expensive searches - that are fundamentally mathematically not guaranteed 
to have converged - we have relied on the negative volume of the Wigner function for our analysis of non-zero temperature 
environments for number states and their finite superpositions. We close with a short conclusion and discussion of our findings 
in Section V. 
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II. THE MASTER EQUATION AND SOLUTION 

Consider a quantum harmonic oscillator linearly coupled to a thermal bath in the interaction picture. Assuming that the bath is 
composed of a continuum of oscillators, that the back-action of the system on the bath is negligible, and modeling the interaction 
in the Markovian approximation, the time evolution of the density matrix p is governed by the master equation IE2I1 

p = l(NL[a?}p+(N + l)L[a]p). (1) 

Here the dot represents the time derivative, the Lindblad superoperator is defined as L[0]p = 20pO^ — O^Op — pO^O, 
and a are the raising and lowering operators for the harmonic oscillator, respectively, 7 represents the degree of coupling of the 
oscillator to its environment, and N corresponds to the mean number of thermal photons in the bath. 

We study this problem in the number basis, where the generally complex numbers (n\p\m) = C n ^ m are the elements of p. In 
this representation equation (Q3 becomes 



]_ ■ 7ii 1 nt 

-C n , m = Vn+ Wm + 1(N + l)C„+i. m+ i - (nN + mN + N + — — )C n , m + y/riy/mNCn-i^-!, (2) 
7 2 

making it clear that the diagonals decouple; that is, C n , m does not depend on C n +i,i+i or Cj+i. m +i, where I, m, and n are 
integers, I 7^ m and j 7^ n M22I1 . This fact allows us to make the substitution m = n + k, k > into Eq. (O, yielding 

1 



-Cn,™+fe = Vn + TVn + k + 1(N + l)C n+1 , n+k+1 

1 



- ((2n + k + 1)N + ——)C n , n+k + y/nVn + kNC n - llV+k -u (3) 

which allows us to map each diagonal to a column vector C k , where k denotes the distance from the principal diagonal. This 
explicit decoupling reduces the dimensionality of the problem from 2 to 1, making it easier to solve. Specifically, we can now 
write the problem as 

-±-C k = A[k]C k , (4) 
7 at 

which is much more tractable than Eq. (O, since it involves the derivative of a vector, rather than the derivative of a matrix. 

A[k] is an infinite matrix and is formally difficult to diagonalize. In Appendix A we show how to use numerical techniques 
to formulate an ansatz (that we go on to prove) for the diagonalization. Solving the problem in diagonal space, and then 
transforming back to the original number space, the general solution is seen to be 



A- 



where dj represents infinitely many constants to be determined by the initial conditions. By using an infinite set of orthonormal 
initial conditions (again see Appendix A) we generate an orthonormal solution basis, 



n n n 



i=0 u=0 

,(fc+2) 7 t/2 mi ^») I !• ( V A ) / N , 1 U 



y ' \3J \U (»+±e*\ (6) 
((N + l)eT* — jyy+v+k+i ^ n + k^ V N J ' 

where m > n. The techniques outlined above and shown in Appendix A can be applied to the case where n > m by making the 
substitution n = m + k into Eq. (O, yielding the almost identical orthonormal basis 
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Equations (|6]l and (0 together allow for a simple construction of a general solution to an initial density matrix p[Q] explicitly 
in terms of its elements represented in number space (whereas using only Eq. (0, one would have to derive the constants aj 
each time). That construction is necessarily piecewise, since the cases n > m and n < to must be treated separately. The final 
solution is 

[0]C n , m [z, t] n < TO 

For a full derivation of this solution see Appendix A. The infinite sums of Eq. ^ converge (assuming the initial density operator 
represents a quantum mechanically valid state) so that any desired accuracy can be achieved by the inclusion of sufficient terms. 

Equation ^ simplifies greatly under certain special conditions of interest. For example, if the initial state is a number state, 
we substitute k = into Eq. (O, and the solution takes the form 



P n [t,n ] = 



Gv + l) 1)Z+J (")(^)(e^(iV + l)-^) i+m S(0(0( N ^ ■ (9) 

where the occupation probabilities P n = C n<n are the diagonal entries of the density matrix. Janszky et. al. derived an equivalent 
result for these initial conditions in [13] by solving Eq. (fl} in the Wigner representation. 

A particularly important simplification we will consider is the case of number states decaying in a zero-temperature bath. In 
this case the solution simplifies greatly: 

P "W = E ( to ) ( n ) (-l)" +m e^ mt , (10) 

m—n ^ / \ / 

which is in agreement with Lu lfl4ll . who pioneered the idea of changing to the number basis, but did so only for the zero- 
temperature case. 

Serafini et. al. derived a solution to a more general form of Eq. (fl~|i which allows for the bath to be squeezed, or phase- 
sensitive lfl5ll . Their solution, given an initial characteristic function x(x,p) i s 

X (x,p) = x{xe-^/ 2 ,pe-^' 2 )e- {xp)a ^ )( 1 "^ 7t )/ 2 (11) 

where is the covariance matrix of the most general single mode Gaussian reservoir. Although formally able to solve the same 
class of equations as our solution, it is numerically simpler to compute the dynamics of measures of non-classicality using our 
number-space representation, rather than working in characteristic space and calculating the number-space entries of the density 
matrix from the characteristic function at each time step. 

Given that our solution is already represented in number space, it is well-suited to comparison of state tomography experi- 
ments. One such experiment was carried out by Brane et. al. lllql . who used rubidium atoms in Rydberg states excited to the 
principle quantum number 50 and 51, and a pulsed microwave source to excite Fock states in a high-Q superconducting cavity. 
They studied the behavior of the Fock states as they decohered in a bath with mean-photon number N = 0.06 using QND 
(quantum non-demolition) measurements. We found that our analytical solution agreed extremely well with their results that 
present the population of the various number states as a function of time. 
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HI. MEASURES OF NON-CLASSICALITY 

In order to quantify how non-classical a given state is, it is necessary to first define what constitutes a classical state. One 
definition of a classical state is one that exists as a possible solution in a classical dynamical system (even if the state satisfies 
the laws of quantum mechanics as well). The formal answers to this question start from this intuition, with the work of Klauder, 
Glauber and Sudarshan for example in considering the coherent states in quantum optic s i23U25ll : it was later shown ll26ll that 
coherent states could be thought of as classical states. The formal criterion provided for a state to be classical was that its 
Glauber-Sudarshan P-representation be positive definite and no more singular than a delta function. 

Another way of determining classical states comes from the intuition developed via the work on classical-quantum correspon- 
dence, particularly as it pertains to classically chaotic systems (see for example e.g. early definitive work by Berry, and later 
work represented by that of Wilkie et. al. 11271 - 13 010 . Berry pioneered the understanding that quantum-classical correspondence 
was best analyzed in the the Wigner-Weyl representation, which he was able to prove is the only phase-space formulation of 
quantum mechanics that approaches the classical limit of Liouville phase-space mechanics. In this case, the criteria for a state 
to be considered classical are that its Wigner representation be continuous, smooth, positive definite, and normalized. Coher- 
ent states are gaussians in the Wigner representation that minimize the uncertainty principle, and satisfy all the properties of a 
classical probability distribution. They are indistinguishable from a classical probability distribution that starts with a minimum 
uncertainty. 

A physically-argued justification that coherent states are the quantum analogues of points in phase-space also comes from 
Zurek et. al. Bill , who showed that coherent states are the most stable in a thermal bath. Specifically they studied decohering 
pure states weakly coupled to a thermal bath, and showed that the coherent states produced the least entropy as they decayed. It 
has also been established elsewhere that the coherent states are the only pure states that are classical [32]. In phase-space this 
argument takes the form of a theorem that the only positive-definite pure state Wigner functions correspond to Gaussian states. 

The most common mixed states which fit all the above descriptions of classical states are thermal states - whether centered at 
the origin or displaced - and have come to be a standard addition to the set of classical states when studying relative measures 
(for example see J^HHHt]). Thermal states centered at the origin are the steady state solutions to Eq. ([T), are defined as 

1 ~ / N \« 

*» = W+i^{n + i) ' (12) 

n— ' 

and represent a thermal distribution of number states, with N defined, as in Eq. (fTJ, to be the mean number of thermal photons 
in the bath. The displaced thermal states are defined using the displacement operator D(a) = expoTa^ — a* a: 

b{a) Pth b- l {a). (13) 

However, as we will discuss in Section IV, including the displaced thermal states in the set of classical states is not sufficient for 
studying the dynamics of the decohering harmonic oscillator, and we will present a novel addition to the basis of classical states. 

We note several caveats here: (a) First, although the quantum optics and the phase-space perspectives agree on the classicality 
of thermal states, there is an important way in which they disagree. Specifically, the phase-space perspective considers all 
gaussians to be classical even though this includes the squeezed coherent states, since the criteria that define a valid classical 
Liouville distribution are invariant under squeezing. However, squeezed states are considered non-classical in quantum optics 
since such states exhibit sub-Poissonian photon statistics. Our perspective is based in phase-space and we do not consider 
squeezed states as a result. 

(b) Of course these classical states - the coherent states and the displaced thermal states are indeed quantum mechanical states 
in that they satisfy the appropriate quantum equations and evolve correctly in the quantum systems. However, their properties 
are indistinguishable from those of classical Liouville probability distributions with a requirement of minimum uncertainty. 

(c) It is also important here to stress that we are only considering unipartite systems. If for example two coherent states are 
entangled (creating a cat state), or a coherent state is used to entangle two other particles, that coherent state can no longer be 
considered as a single entity, and the total multipartite system would certainly not be considered classical. 

Having established a set of classical states, it is then possible to classify how non-classical another state is by considering 
how different the given state is from the set of classical states. We will consider three of the standard ways of quantifying this 
distance between an unknown state and the set of classical states. 

These are relative measures, comparing unknown states to those known to be classical. These are formally well-defined. 
However, as we see below, the non-classicality of a state is defined in terms of the infimum (or in the case of classicality 
measures, the supremum) of the distance from all classical states. Given that it is impossible to search through this infinite 
set of all classical states, there is a certain amount of ambiguity about these relative measures. Progress can be made using an 
incomplete basis, however, as has been previously achieved, and in particular in some cases upper and lower bounds on the 
non-classical distance can be established yj]. 

Other ways of quantifying non-classicality can be understood as absolute measures in that they depend only on inherent 
properties of the state under consideration, and so avoid the problem of defining a set of classical states. Perhaps the simplest 
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and most well-known examples of this type of measure is the purity, Trp 2 , or its close relative, the 2-entropy, Log[Trp 2 ]. For 
example, as mentioned in passing above, Zurek et. al. i3lll studied dynamics by calculating the entropy (specifically the 2- 
entropy) generated by the decohering pure states, and offered further evidence that the coherent states are the most classical pure 
states. Marian et. al. have also used the same measure to show that in most cases as Fock states decay in the Markovian bath, 
they reach a point of maximal mixing, after which time the purity or 2-entropy increases again. This result indicates that while 
the purity can be a reasonable rule of thumb in the static case, it does not hold up in the dynamical case as a sensible measure 
of non-classicality. This is because for an open system it does not conform with intuition that non-classicality should decay 
to a minimum, and then increase again under the effects of decoherence. Further, coherent states, which have been defined as 
the most classical states of all by many different authors, have unit (maximum) purity. Dodonov et. al. make similar arguments 
against using purity as a measure of non-classicality yfl. 

Many absolute measures involve representing the uncharacterized state in phase space and measuring some property of that 
representation. Paavola et. al. , who studied dynamical systems using nonclassical measures, present five such measures and 
use them to study the decoherence of cat states 111 211 . They consider how long it takes for the Wigner function to become fully 
positive, the time it takes for the Glauber-Sudarshan function to become fully positive, the time it takes for the interference fringe 
in the Wigner representation to disappear, the Vogel criterion, and the Klyshko criterion. They show that of these five measures, 
all but the interference fringe technique show that the cat state becomes classical in a finite amount of time, and that those times 
are all on the same order of magnitude. Since they have already shown these four approaches to give similar results, we will 
only consider one of the measures in this paper. The only other paper studying dynamics of non-classicality that we are aware of 
also studied the negative volume of the Wigner function of cat states [ 34] . In particular, that paper presented analytical formulas 
for when the Wigner functions of cat states become fully positive. 

With this as background, the four measures we present results on below are 

1. The first is a measure discussed in detail by Hillery, being defined as the infimum of the trace norm of the difference 
between the quantum state and the set of all classical states p c i [1]: 

t]h[p}=^\\p- Pc\\ (14) 

Pal 

where P| I = Tr[VX4J] is the trace norm of a matrix. Hillery used this measure to prove that there will always be a finite 
distance between non-classical states and classical states. He also used it to study the harmonic oscillator eigenstates, and 
proved that t\h is bounded from above and below by 

1 - In <VH< 2y/l - 7 n (15) 

and that j n must decrease at least as quickly as iiT 1 !" 1 for large n. This result simply rests on the existence of a classical 
basis, and does not require its construction. In Section IV. we use a partial classical basis to show that non-classicality 
in fact increases monotonically and asymptotically to the maximum value predicted of 2. Though it is not immediately 
evident from Eq. (fl4l >. the maximum value of t]h is in fact 2, meaning that increasing number states approach the maximal 
value of non-classicality. 

2. The second measure, the Bures distance, uses the Uhlmann fidelity, given by B[pi, P2] — Tr[^/ \/pip2^/pi}, which was 
shown to be the density operator generalization of KV'ilV^)! by Uhlmann Q3|] _ The actual Bures distance is calculated 
using Jl] 

tjb = inf y/2(l-\B\p, Pcl ]\). (16) 

Pal 

This measure was used to determine analytically the non-classicality of squeezed, displaced thermal states [2], and later 
the non-classicality of two such states when entangled l33ll . In Section IV. we show that this measure agrees well with the 
Hillery measure in the case of the number states. The only real difference, that t]b increases asymptotically to \/2 instead 
of 2, is due simply to its construction. Eq. ( fTBI l clearly shows that the maximum value rjs can attain is \[2. 

3. The third measure (the Dodonov measure), is a measure of classicality, rather than non-classicality, as it computes how 
close, rather than far, a state is to the classical basis. It is presented in [3], and is a modified form of the overlap of density 
operators: 

T)n[p] = supTi[p'p' d ] (17) 

Pal 

where p' = . p is the density operator renormalized by its purity. The purpose of this renormalization is to make 

y/Tr[p 2 ] 

it easier to compare mixed states. Notice that for mixed states Trf//^] = 1 if Pi = P2, whereas Tr[pip2] < 1 in the 
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same case. Thus, the renormalized overlap allows one to easily see how similar two density operators are, regardless of 
the degree of mixedness. Dodonov et. al. used this measure to show that the classicality of the number states decreases 
monotonically and asymptotically to 0. Again this means that increasing number states approach the maximal value of 
non-classicality. Specifically they showed that rjn = er n n n fn\, which, using the Sterling approximation, reduces to 
r)r> s» (27m) _1 for large n. 

4. The measure we will consider from the family of absolute measures is the negative volume of the Wigner function, or 
negativity, 

/OO /"OO 
/ (\W[x,p}\-W[x )P })dxd P> (18) 
-oo J — oo 

which is particularly appealing because it is easily calculated using numerical integration. It is also consistent with the 
intuition resulting from the study of quantum-classical correspondence in phase-space following Berry. Kenfack et. al. 
have already used r\w to study the harmonic oscillator eigenstates (for the closed, static problem) Their results agreed 
very well with those of Dodonov et. al. in that Kenfack et. al. found that 77^/ increases like y/n fl. 

We now present results using these measures, starting with the previously well-studied (see e.g. Ollltl) base-line case where 
there is no coupling with the environment (which can alternatively be understood as the initial condition for all the cases). 



IV. RESULTS 



A. No Environmental Coupling 

As discussed in section III, the Hillery and Dodonov et. al. measures have been previously used to show that the non- 
classicality of Fock states increases with increasing n in the absence of any environmental coupling. Dodonov et. al. derived 
an analytical expression for this increase using the classical basis composed of the displaced thermal states. In studying this 
problem we found and present here an improvement upon both these basis, which minimizes r\u further (recall that this measure 
seeks the infimum over all classical states), and further maximizes t)h and r\B (recall that these measures seek the supremum 
over all classical states). 

We discovered this new basis - particularly of relevance to decohering states - through our own numerical experiments as 
well as from studying previous results on the behavior of initial number states [f35i] . What has been observed is as follows: the 
Wigner functions for the number states start as Laguerre polynomials in phase-space, with multiple non-classical fringes. As the 
states evolve in the presence of the environment, they lose their initial fringes until they become slightly broadened versions of 
the classical microcanonical states. That is, they become positive-definite distributions which are sharply peaked in energy space 
on the appropriate classical energy, and are otherwise evenly distributed along the classical orbit. The classical microcanonical 
state at a given energy E is, of course, 5(E — H(p, q)); the states we are talking about are not quite as singular as these delta 
functions in energy in having a slight spread in energy. 

These 'thermally-broadened' microcanonical states are a natural choice as classical states which might be closest to these 
quantum eigenstates. In studying quantum-classical correspondence averaging over neighboring states has been shown to be 
necessary, yielding states similar to these thermally broadened micro-canonical states 027142911 . Irrespective of the intuition, as 
we show below, empirically these states work very well indeed. 

We represent these states by 

00 

p+ = Y / an P +, (19) 

where X^^Lo a « = 1 ar >d a n <= [0, 1], and p+ represents a number state that has evolved according to Eq. (Q]) to a time 
t* = - log ^34lh wn i cn guarantees that it is positive definite in the Wigner representation. This is an extremely large set 

and is impossible to explore completely. For our purposes, we found that the most useful subset of p + is 

p+ = (x + 1 - v)pt + {v- x)p+ +l (20) 

where x = \nt\y\ is the truncated integer of v, and v varies continuously. For example by this definition = .\p^ + -9p^. 

In Fig.Q]we show the results of computing these measures of non-classicality for eigen-states of the harmonic oscillator using 
multiple measures, and using three different classical bases: the set of all coherent states, the set of all thermal states, and the set 
p+. We note that although clearly p+ outperforms the other bases, it does not provide any new insight into the behavior of the 
system. 
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FIG. 1. For the static (or closed/non-interacting) case we show different measures of non-classicality for eigenstates of the harmonic oscillator. 
Specifically we show the (a) Hillery Distance, ijh, (b) Bures Distance, ijb, and (c) Dodonov Overlap, 770, each computed for three different 
classical bases, and (d) negativity, r\w ■ 



lH IB 




(c) (d) 



B. Zero Temperature 

The results of the previous section seem to violate the correspondence principle on the face of it, since they show higher n 
states as more non-classical in contradiction to the rule of thumb that higher number states should be less non-classical. This 
apparent paradox has in fact been well known for decades: that higher n states display more non-classicality - via rapidly 
oscillating fringes, for example - in violation of the correspondence principle intuition. It has been previously shown that 
the correct way to approach this issue is by averaging the Wigner function over a small energy spread or phase-space spread 
(usually evoked as a resolution limit). When such averaging is done for the Fock states it appropriately smooths the fringes, 
or oscillations, of the number states 11271 12811 . The same smoothing effect occurs explicitly when studying dynamics. As soon 
as environmental coupling is considered, the number states evolve over time into states resembling the micro-canonical states 
(in this case the 'averaging' as being done by the environment). Higher n states are much more sensitive to decoherence than 
lower n states, which means smoothing occurs much more rapidly for higher n states, and yields the appropriate correspondence 
principle intuition. 

Another way of understanding this idea, that the correspondence principle requires averaging (either mathematical or envi- 
ronmental), starts from the notion that the behavior of closed systems is singular and non-physical, in that states in such systems 
evolve completely independent of their environment; they do not even interact with the vacuum. This behavior does not survive 
when even the smallest amount of interaction with the environment is introduced. 

Specifically, for example, by introducing a zero-temperature bath into the problem, the behavior of the system changes 
markedly from that displayed in Fig. Q] Using the expanded classical basis outlined above, each of the three relative mea- 
sures, along with the absolute negativity measure, yield results in agreement with the correspondence principle (see Fig.|2]l. That 
is, all of the non-classicality curves shown for finite jt peak, and then asymptotically decrease to with increasing n (although 
this is harder to see for the Hillery and Bures measures. The new basis turns out to be critical in detecting this. Without it (that 
is, using only the displaced thermal staes), according to the relative measures the non-classicality of the decohering eigenstates 
shows a monotonic increase for all yt curves. 

Though each of the measures show the same general result, each one differs somewhat on the details of the behavior. Ac- 
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cording to (c) the Dodonov overlap and (d) negativity, this non-monotonicity is relatively straightforward. Each state decays 
smoothly, and the curves for each state each other state only once, so that the non-classicality peak decreases in n monotonically. 
However, the (a) Hillery and (b)Bures distances tell a more complicated story. As can be seen in Fig. |3(a)| individual eigenstates 
do not decay smoothly as a function of time according to the Hillery distance; specifically the 4th and 7th excited states exhibit 
multiple corners, giving rise to multiple crossovers. This naturally leads to the more complex behavior seen in Fig. |2(a)| The 
time curves in this case display not only a maximum value, but also local minima. 

The Bures distance shows individual eigenstates decaying more smoothly than the Hillery distance, though Fig. |2(b)| does 
show occasional local minima. The more unusual signature present in this measure is that many of the eigenstates never cross 
the decay curve of the 1st eigenstate, as seen in Fig. |3(b)| For this reason the non-classicality peak never reaches n = 1. The 
lowest value it reaches n — 2, after which time it begins moving up again. 

These unexpected features may be explainable by the fact that we used a classical basis that is far from complete. Ideally the 
Hilbert-space measures rely on searching over the entire set of classical states, which is practically impossible. Unfortunately 
this inability to represent the complete classical basis can lead to incorrect results. Indeed, as we have mentioned, when using 
the coherent or thermal states as reference states the non-monotonic signature disappears altogether. It may be reasonable to 
assume that the non-classicality peak would come to rest at n = 1 according to the Bures distance as well if a more complete 
basis were used. Similarly, it is arguable that the Hillery distance would decay smoothly for all harmonic oscillator eigenstates 
if a more complete basis were found. 

Independent of these speculative arguments about these unusual features, all four measures show qualitatively similar behavior. 
That is, in all cases, an initially ('yt = 0) monotonic dependence on quantum number transforms into non-monotonic behavior 
in the presence of decoherence. 

We comment here about the relative computational difficulty of the different measures. We have just demonstrated that in 
the case of zero temperature negativity (the negative volume of the Wigner function) shows good agreement with the other 
three measures discussed. However, comparing these measures in more general cases is even more formidable challenging. 
Specifically we would have to search even more broadly to discover the best classical bases for each problem. Further, the 
zero temperature problem is relatively convenient computationally speaking, because it deals with finite density matrices in the 
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FIG. 3. The decrease in non-classicality as a function of time for the 1st, 4th and 7th Fock states as measured by the (a) Hillery distance, (b) 
Bures distance, (c) Dodonov overlap, and (d) negativity. 




number representation. The non-zero temperature leads to infinite density matrices in number space, and makes calculating non- 
classicality more and more expensive as temperature increases. For negativity, only one numerical integral must be calculated 
at every time step, whereas for the other measures hundreds of values must be calculated and compared at each time step, so 
the difficulty of the problem scales much more quickly. Given the computational difficulty of the relative measures, and the 
uncertainty of finding adequate bases for further study, we choose to use solely the negativity as a measure of non-classicality 
for the rest of our studies. 

For those who are still suspicious of using negativity as a measure of non-classicality, these results may well seem unsatisfying. 
However, we would like to remind the reader that the only other two studies of dynamics we found lfl2[ l34ll dealt solely with 
phase-space measures. We have attempted to validate this approach by comparing the measures where we found it possible, 
just as Kenfack et. al. did before [4] by considering the static case. In addition, we point out a further pragmatic argument for 
considering negativity to be useful in quantifying non-classicality, which comes from somewhat recent discoveries in the field 
of quantum computation that deal with a discrete analogue of the Wigner function. Cormick et. al. have shown that the only 
non-negative pure states in this representation are the stabilizer states [36]. Galvao has demonstrated the same result, and used it 
to support the conjecture that negativity of the discrete Wigner function is a necessary condition for the exponential speedup of 
quantum computation with pure states IB7I1 . 
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C. Non-zero Temperature 



In Fig. |4(a)| we show negativity versus Fock state for different values of jt at finite temperature (N = 0.06). Notice that at 
7< = the negativity increases monotonically in n, but that at any nonzero time there is a peak in negativity, again confirming 
our argument that the monotonicity of negativity at = is singular and unphysical. 

This is reinforced by Fig. |4(b)| which shows evolution of the lowest 16 harmonic oscillator eigenstates to -ft = 0.15 for 
different bath temperatures. Here, as we already know from Fig. |2(d)| even the N — curve has a peak, in contrast with the 
t = peak in Figs. |2(d)] and |4(a)| We see that temperature has a somewhat similar effect on non-classicality as time in that higher 
temperature corresponds to greater decoherence, just as greater values of jt correspond to greater coherence. On the other hand, 
the mere fact of t > leads to a peak, and the value of jt places an upper limit on what value of n the non-classicality peak 
takes. 



FIG. 4. (a) The time dependence (at constant bath temperature) of the negativity of the Wigner function for Fock states evolving in a bath of 
mean photon number N = 0.06. Each curve is a snapshot at different values of time of negativity as a function of eigennumber n. (b) The 
temperature dependence (at constant time) of negativity of the Wigner function as a function of n. Each curve is a snapshot at different value 
of the bath temperature, of the negativity as a function of eigennumber n, all taken at the same unitless time jt = .15. 




(a) 




^— N=0 

N=0.06 
^— N=0.25 

N=0.5 
-^N=l 



n 
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D. Non-diagonal Initial Density Matrices 



As a final, brief foray into the rich landscape of potential states that may be explored, we considered four other types of 
initial conditions (well-suited to the form of our solution), all evolving in a Markovian bath of mean-photon number N = 0.06, 
which, as shown by Brune et. al. is an experimentally realizable temperature II 1611 . The first, displayed in Fig. |5(a)| are of 
the form ip n [t = 0] = -j= (\n — 1) + |n)), and represent the simplest form of qubit states. The second, displayed in Fig. 

|5(b)| are of the form ip n [t = 0] = A= (\n — 1) + \n + 1)), and represent less conventional qubit states. The third, displayed 

in Fig. |5(c)| are of the form ip n [t = 0] = -4= Ylm=l \ m )- The f° ur thi displayed in Fig. |5(d)| are of the form ip n [t — 0] = 



En — 1 
m— 1 



2^1 



1 2m- 1) 



2("-!)/ 2 



\2n — 1). This last two are brief forays into the large expanse of possible qudits. 



FIG. 5. These plots show "non-classicality" as a function of initial condition n in a bath of mean-photon number N = 0.06 at specific values 

of ft for states described by (a) ip n [t = 0] = -j= (\n — I) + \n}), 

(b)V„[* = 0] = J=(|n-l) + |n+l»,(c)V„[t = 0] = ^£^ =1 |m),and 

(d) i> n [t = 0] = E™=\ ^|2m - 1) + ^m\1n - 1). 




f=f^»[0] o.o 



f«M0] o.o 




These results clearly show that in general "non-classicality" is not as straight-forward for general states as for Fock states 
might at first suggest. However, there are elements of similarity to the Fock states. Particularly in the case of Figs. |5(a)| and 
|5(b)| we see a monotonic rise in "non-classicality" at t = and a "non-classicality" peak at all other times (until all states 
have transitioned to classical). The major difference lies in the fact that the -^|0) -I- ^75 1 1) state never appears to cross the 

^7||1) + 1 2) state; nor does the ^|0) + ^= j 2) state appear to cross the + ^ |3) state, so that the "non-classicality" 

peak appears to stop decreasing at n = 2. This is not entirely surprising, since each of these n = l states are constructed partially 
from the ground state, which is by our definition a classical state. 

The non-monotonicity of the t = curve is somewhat unexpected in Fig. |5(c)| However, the time-dependent behavior is fairly 
similar to the Fock state case, in that ultimately the "non-classicality" peak occurs at n = 1 (the state corresponding to the initial 
condition \tp) = |1)). In contrast, the non-monotonicity of the t = curve in Fig. |5(d)| is not so surprising. It makes sense that at 
first as the states become more complex, the "non-classicality" increases. However , ultimately these states are approaching an 
infinite superposition of eigenstates, rather similar to a coherent state, and therefore would be expected to approach a classical 
state. Here again the "non-classicality" peak behaves somewhat unexpectedly. Though it eventually comes to rest at n = 1, it 
skips n — 2 entirely. 

These last two may have important implications for understanding what types of states will be useful for qudit encoding, 
especially since negativity is believed to be a necessary condition for quantum computation speedups Il37ll although it is the 
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negativity of multiple entangled qudits that is specifically important. However, it may be useful to understand the behavior of a 
single qudit first. 



V. CONCLUSIONS 

We have presented above a general solution to the harmonic oscillator linearly coupled to a Markovian bath, and have studied 
the behavior of four different measures of non-classicality using the Fock states and special finite superposition of Fock states. 
By introducing a new classical basis, we were able to use the three relative measures to show that this system obeys the corre- 
spondence principle and indicates a non-monotonic transition from quantum to classical. We also showed that negativity gives 
the same result. 

The result obtained using the three relative measures relies upon the new classical basis we introduce: Without this new basis, 
the measures appear to contradict the correspondence principle. We find that negativity is arguably a more versatile and useful 
measure for studying the dynamics of systems - particularly where an appropriate classical basis is harder to guess. Finally, we 
used negativity to explore the dynamics of number states decaying in a non-zero temperature bath, and then generalized to a few 
finite superpositions of eigenstates; the last serves as a brief survey of the rich landscape of states that are still open for study. 
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Appendix A: Deriving the Solution 



We solved the master equation for the simple harmonic oscillator linearly coupled to a Markovian bath of oscillators, given in 
Eq. [TJ by representing that equation in the number basis (see Eq. O, and then making the substitution m = n + k, which leads 
to Eq.[3] and is reproduced below: 



1 ■ 



C n ,n+k = \/n+lVn + k + l(N + 1)C„ +1 , n +fc+i 



((2n + k + 1)N + 2n ^ k )C n , n+k + VnVn + kNCn^n+k-L 



(Al) 



As previously stated, Eq. lAll is equivalent to the much more compact vector equation = A[k]C k , where A[k] is an infinite 
2-dimensional matrix and C k is a vector composed of the elements C n . n +k- Since A[k] is an infinite dimensional tri-diagonal 
matrix, we cannot diagonalize it analytically. Instead we search for the eigenvalues numerically by finding the eigenvalues of 
truncated matrices A[k]d of dimension (d + 1) x (d + 1). By graphing the eigenvalues as d becomes large, we see a very 
compelling convergence to Xj = — (j + |), where j denotes the jth eigenvalue and spans the nonnegative integers. 

We then take this premise as our Ansatz, and proceeded to find the eigenvectors using these eigenvalues. We compute a finite 
number of entries for multiple eigenvectors until a pattern emerges. We see that the Ith entry of the jth eigenvector is 




(-1)' 



2 = 



N~ 



1 



k + i 
k 



(A2) 



which we use as our second Ansatz (note that we have as yet not proven the first). 

We then construct finite truncations of the transformation matrix T from Eq. (IA21 > and invert them to search for the analytical 
form of T _1 . Again a pattern emerges, suggesting that the Ith, jth entry of T _1 is 




N - y>i r 

i/=0 



1 



k + v 
k 



(A3) 



We prove that this is the inverse transformation matrix in Appendix B by showing that Eq. dA2b multiplied by Eq. ( IA3b yields 
the Kronecker delta. 

Using the eigenvalues from our first Ansatz, we construct the diagonal matrix D = T~ 1 A[k]T, which is of course comprised 
of the eigenvalues shown along the diagonal. In diagonal space, we have the differential equation ^^affc = Dx, where Xk 

TCk- The solution to this set of differential equations is trivial: 



(A4) 



and can be expressed in vector notation as x k = e lkt l 2 Ea, where 

/. 



E 



-njt q 



e" 2 ^ 
e~7* 
1 ) 



(A5) 
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and a is composed of the unknown constants a n . 

Transforming back into the original space, we have Ck [t] = Te~ lkt / 2 Ea, or 



<w - X>y ( "X " ) G£r)" £<-»* ( I ) ( J ) <"> 

The constants a n depend upon the initial conditions, and can be solved by setting t = 0: Ck[0] — Ta, or a — T~ 1 Ck[0)- We 
now have all the tools we need to write down the final solution in vector notation: 

C k [t] = fe-^ 2 Ef- 1 C k [0] . (A7) 

In order to find the closed form expression for an entry of a general density matrix C n>n +k< we define the initial conditions to 

be 

Cn, n+k [0) = {l] n n ^2 , (A8) 

which gives rise to Eq. © upon substituting them into Eq. (IA71 >. Since these initial condition vectors represent an orthonormal 
set, Eq. (O can be used to find the solution to any initial condition. This idea is represented mathematically in Eq. (0. 

However, this solution rests on two Ansatze: that the eigenvalues of A[k] are Xj — —(J + |) and that the eigenvectors are 
given by Eq. (IA2I) . These Ansatze are proven by substituting Eq. (0 into Eq. (f3]), and showing that the solution we found does 
in fact solve the original differential equation (see Appendix C). Since there are sufficient a n to span the entire solution space, 
the uniqueness and existence theorem guarantees that our solution is correct, and by extension that our Ansatze were correct. 

Appendix B: Inversion Proof 

In order to prove that the matrix, here denoted B, given by equation (I A3t and reproduced below 



(N , ... 

V ' !/=0 



is the inverse of T, given by equation (IA2t and also reproduced below 



^?- ir (;) r ife (bi) 



k + j\fk + l\ N l ^, ..„ / l\ AT _ V ( j\ 1 



k 



k 



we will show that the ath row of T multiplied by the bth column of B is the Kronecker delta: 



k + a\fk + b\ N a \^v^, , -, T -i-v ( a \ ( b \ 1 



, Wf-iv +i/ jv^, , , , 

k J (N + 4-i ^ \i \ v { k + v\ { k + i 



k \ k 



E 



m \ I m \ I m 



k\ ( N 



(B3) 



i J \ v J \ k J \N+l 
Dealing with the infinite sum in square brackets (letting (^j^fjj = r), it is easy to show that 

E (?) (?) ( TO fe fe ) rm = feT^ rfc {T!^((l-V 1 )}" (B4) 
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We first simplify the term in curly braces. Using the fact that: 



v=0 

and setting f[r] = (1 — r) - ' -1 and g[r] = r l we have: 

i 



|r(/H5M) = E ( i ) f {l ~ v) V\9 {v) \r], (B5) 



«dr» V(l-r)*+V "^V"/ *l (*-01 ' 

which simplifies to 



(1 - rV'+ i+1 ^ V f 



By the binomial expansion this becomes: 



Reordering the sum and grouping by r M we have: 



r l+l ^ -u,^ ( I \ ( i + l ~ v \ ( v 



The coefficients of can be simplified by summing over v. We first reindex the sum by setting v = fi + k: 

^(^('^-^(^(-ir— g( M it)(' +, 7"*)(' , ;*)(- 1 ) M *- < Blm 

The right-hand sum can be rewritten as a multiple of a hypergeometric series: 

» ( 't|l~^;, r 2 J\[-a-M),M-*;M-»-*;i], (BID 

allowing us to use the Chu-Vandermonde identity, which states that 

2 F 1 [-n,/3; 7 ;l]= (7 ~^" ' < B12 > 

(7)n 

where (x)„ = x(a; + l)...(x + n — 1); n > is the Pochanner symbol. Substituting this solution for the hypergeometric series 
yields 



v-f I \fi + l-n-k \(fi + k\ r _ lV ,+„ = (i + (-Ot-M 



which can be rewritten as 



(i - n)\(i - fulfil n\(i + i-n)\ (i- - • 

Therefore 



(B13) 



(B14) 



EC)( "',-")( ;<-')" + "=: HI) <B> 5 , 
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and 



m=0 



X ( "! ) ( " ' I r 



r l+l J^f I \ l i 



(l-r) 



i/=0 



Returning now to Eq. SB4\ . we have 



1 d fc 
fc! dr fc 



„i+i> 



tE 



(1 - r) i +^+ 1 ^ \ 3 J \ 3 

3=0 K ' v 



which can be written as: 



fc! ^ V 7 V i 7 V M 7 rfr fc -^ 



d' 1 



This becomes 



r \— s I I \ I V 



fe+1 E 



fl — r-) i + l/ + fe + 1 ^ ' V 7 / V 7 



E 

Ai=0 



i + ^ + fc-/i^ ( i + v + k — j\ (l-r 
i + v 



and by the binomial expansion we have 

y - V+k+l E 



Ji+v+k v 



(1 - r) 1 



3 J \ 3 



r - 3 ^ i + v + k — /z ^ ^ i+v+k-j^ ^ M 



3=0 x ' x ' n=0 ^ ' v 

By rearranging the sums over w and 7X this expression can be rewritten as 



tE 



n — r)«+H-fc+i \ 3 J \ 3 



E 



E 



i + ^ + fc-7i^ f i + i/ + fc - j ^ ^ /1 ( 



i + v 



The term from (1B2U in square brackets is actually a hypergeo metric series, and can be written as 

(i + v + k- + v + k - w)\ . . . 

-7— rp- rrr— TTT 2-Fl ~ fc — U>), (u> + J — fc — Z — Z^); — fc — Z — ^); 1 , 

(« + ^)!w!(fc — + u> — fc)! 
which allows us to apply the Chu-Vandermonde identity again to arrive at 

- ( 1 + 

i + v 



E 



/i 



i+v+k- ^\ ( i + v + fc - j \ ( fi \ r^xft+w _ ( i + v + k- j 

w I ^ ' \ w j \ k — w 
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(B16) 



(B17) 



(B18) 



(B19) 



(B20) 



(B21) 



(B22) 



(B23) 



Therefore 
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m=0 



m\ ( m \ ( m+k 



(B24) 
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1+1 E 



3=0 



E 
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i + v + k — j 
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k — w 



Recalling that r = j^xj-, equation ( IB3t can be written without the infinite sum as 

a+k a b u k 
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k 



N 
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EEEEM) i+ i: 
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(B25) 
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N 



3 J \ 3 J \N + 1 
Rearranging the sums over j and w yields: 



i + v + k - j \ ( j \ ( N 
w \k-w \N + 1 



k J \N + 



\ a+k a b u 
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The term in square brackets in Eq. (1B271 is a 3^2 hypergeometric series: 
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to=0 
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(B27) 



3 J \ 3 
which can be simplified to 
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by use of the Pfaff-Saalschutz identity, which states that 



3 F 2 [—71, a, b; c, 1 + a + b - c — n; 1] 



(c - a)n(c - b) r 
(c)„(c - a - 6), 



Thus equation ( IB3b reduces to 
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Rearranging the terms for clarity, we have: 
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Calling the combinatorial identity for the kronecker delta 

E(- 1 ) m+J 

twice we have 
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Therefore T£? = /, and we have proven that B — T 1 . 
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Appendix C: Uniqueness and Existence Proof 



We prove that equation ©, reproduced below 



is the solution to equation ([3), reproduced below 



-C'n.n+fc = V" + iVn + fc + 1(JV + l)C n +l,n+fc+l 

7 



2n + k, 



- ((2n + fc + l)iV + — — )C„ in+fe + yWra + fcJVC n _i in+fe _i 

by substituting equation ( ICU into equation (IC2l i and showing that both sides are equal. 
We start by breaking the righthand side of equation (IC2t into three parts: 



(y n +lVn + k + l(N + l)C n+ i, n+fe+ i - (n + fc + l)NC n>n+i \ 



+ {-n(N + l)C„,n+fe + Vny/n + kNCn-^n+k-i} 
Comparing this to the lefthand side 



2 ^TT^n+fc 



j'=0 V v 7 v 7 i=0 



we see that the term in square brackets in equation (IC3t also occurs on the lefthand side. 
Next we consider the term in parentheses from equation (1C3V 
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Using the binomial identity 



n \ / n — 1 \ f n — 1 
k I ~\ k ) + { fc- 1 



we subtract equation (IC6b from equation (IC5t and reindex the sum by i — > i + 1 to get 
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Now we consider the term from equation (1C31 in curly brackets 
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We subtract Eq. ( !C9t from Eq. (IClOt . once again using Eq. (IC7t and reindexing the sum by i — > i + 1, to get 
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If we subtract the lefthand side of Eq. (lC2t from both sides, the term in Eq. (|C3t in square brackets cancels, and the remaining 
terms are 
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From here we examine the series over i. Recognizing that the constant term is zero, the series can be rearranged to be 
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(C13) 



It is now easy to show that each coefficient of (— j^) 1 is 0. Therefore Eq. (1C12| > reduces to = 0, which is equivalent to 
the statement that both sides of equation ( !C2t are equal. Therefore we have shown that equation ( IC1) is a solution to equation 
(1C21 >. Since equation ( ICU has an undetermined constant for each initial condition C n , n +k [0], this proves by the Uniqueness and 
Existence theorem that we have found the complete general solution to equation ( IC2b . 



